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Abstract 

We investigate actions for dynamically triangulated random surfaces that consist of a gaussian or 
area term plus the modulus of the gaussian curvature and compare their behavior with both gaussian 
plus extrinsic curvature and "Steiner" actions. 



Considerable effort has recently been devoted to exploring modifications of the discretized Polyakov 
partition function Q] for a random surface 



T 



N-l 



Y[dX?exp(-S g ), (1) 



where the sum over triangulations means that we have, in effect, a fluid surface. 5^ is just a simple 
gaussian action 

s s = ~ £(*f-*jT, (2) 

<ij> 

where the X's live at the vertices of the triangulation and the sum < ij > is over all the edges. Earlier 
work ||] had made it clear that this action (and variations, such as area and edge length actions [|[ ||) 
failed to lead to a sensible continuum theory because the string tension did not scale so, inspired by 
analytical work on QCD strings and biological membranes 0, an extrinsic curvature or "stiffness" term 
was added 

S e = (1-Ui- Uj), (3) 

where n^, rij are the normals on neighboring triangles A^, Aj. Simulations of S g + XS e , called the gaussian 
plus extrinsic curvature (GPEC) action, seemed to indicate that there was a second order phase transition 
(the "crumpling transition" ) from a small A crumpled phase to a large A smooth phase at which one might 
hope to define a non-trivial continuum theory [pi. More recent simulations of larger surfaces suggest, 
however, that the transition is not second order |7], ^] [] but that the string tension may possibly still be 
scaling correctly. 

Simpler spin systems on dynamical triangulations, such as Ising and Potts models [fiol , provide some 
reassurance that a non-trivial theory may be lurking at the crumpling transition because they display 
third order transitions and still have a sensible continuum limit. Nonetheless, the rather murky be- 
havior seen in (7[ |s) prompts the question of whether actions with a sharper phase transition can be 
found. One geometrically appealing suggestion was the Steiner action put forward in JHJ, and simulated 
microcanonically in |fl2f 

S st einer = l £ \X? - X?\8( aij ), (4) 

<ij> 

where 0{otij) — \ir — ctij\ and ot^ is the angle between the embedded neighboring triangles with common 
link < ij >. It was pointed out in jl3) that the grand canonical partition function diverged for S steiner 
alone, so we conducted some exploratory simulations of various actions combining edge-length, area or 
gaussian terms with S 'steiner finding particularly sharp transitions for Area + XS 'steiner and S g + XS 'steiner 
actions As this initial work is on the same small triangulations that indicated a second order 

transition for the GPEC action we cannot claim this is strong evidence for a sharp transition with the 
Steiner actions, without further results from larger lattices [ fl5f . 

In this paper we will explore a further possibility for a random surface action which, like the Steiner 
term, is a natural object to consider from a geometrical point of view. For a curve C embedded in three 
dimensions it was shown by Fenchel that 

- / \n\ds > 2 (5) 
n Jc 

where k is the curvature and the equality holds when C is a plane convex curve Jl6[ . For a surface M 
imbedded in three dimensions the Gauss-Bonnet theorem tells us that 
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KdS = X (M) (6) 

M 



where K is now the Gaussian curvature of the surface and x 1S the Euler characteristic. As this is a 
topological invariant it tells us nothing about the configuration of the surface. To get the equivalent of 

1 Unlike the GPEC action on rigid lattices, which almost certainly is [Bl . 
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Fenchel's equ.(||) we take a modulus sign in equ.(||) and find 

\K\dS > 4 - X (M) (7) 
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with the equality holding when the surface is imbedded as a convex surface in three dimensional space. 
This term discretizes to 

S« fl ht = ^|27r-5^^ i | (8) 

' iOO 

where the outer sum is over all the vertices of the triangulation and the inner sum is round the neighbors 
j of a node i. (fiij is the angle subtended by the jth triangle at the ith vertex - see Fig.l. The behavior 
of Sught was measured for the GPEC action in || (rather than being included in the action), where it 
was found to correlate closely with the extrinsic curvature - dropping sharply in value at the crumpling 
transition. This raises the hope that including Sught in the dynamics may be sufficient to give a crumpling 
transition without the assistance of an extrinsic curvature term. 
We will examine the phase structure of both Si = S g + XSught 



Si = \Y.^- X j) 2 + X Y.\^-H^\ ( 9 ) 

and S 2 = Area + XS Ug ht 
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<ij> i 



S 2 =J2^+xJ2\2tt-J2^\ (10) 

A i 

where Aa is the area of a triangle A, in what follows. To do this we employ our by now standard set 
of observables. We included a local factor in the measure for compatibility with our earlier simulations 
which can be exponentiated to give 

S m =^logfe), (11) 

i 

where is the number of neighbors of point i, and d = 3 dimensions. We thus simulated Si, 2 + S m . We 
measured < S m > and the mean maximum number of neighbors < max(qi) > to get some idea of the 
behavior of the intrinsic geometry. The extrinsic geometry was observed by measuring < Sught > and 
its associated specific heat 

C = ^ (< Sf igM >-< S tight > 2 ) (12) 

as well as the gyration radius X2, a measure of the mean size of the surface as seen in the embedding 
space, 

X2 = mhi)Z^- x ^ mj - (13) 



For Si the expectation value of S g can be shown to be d(N — l)/2 by exactly the same argument that is 
used to give this result for the GPEC action, because Su g ht shares the scale invariance of the extrinsic 
curvature term in the GPEC action. This serves as a useful check of equilibration in this case. A further 
useful check is provided by removing the modulus sign in equ.(||) which should then sum to give 47r for 
every configuration by the Gauss-Bonnet theorem. 

The simulation used a Monte Carlo procedure which we have described in some detail elsewhere JlTj . 
It first goes through the mesh moving the X's, carrying out a Metropolis accept/reject at each step, and 
then goes through the mesh again carrying out the "flip" moves on the links, again applying a Metropolis 
accept/reject at each stage. The entire procedure constitutes a sweep. Due to the correlated nature of 
the data, a measurement was taken every tenth sweep and binning techniques were used to analyze the 
errors. We carried out 10K thermalization sweeps followed by 30K measurement sweeps for each data 
point. The acceptance for the X move was monitored and the size of the shift was adjusted to maintain 
an acceptance of around 50 percent. The acceptance for the flip move was also measured, but in this 

2 As a particular case of the theorem: Let M™ be a compact oriented C°° manifold immersed in i? n + JV , such that the 
total absolute curvature equals 2. Then M n belongs to a linear subvariety of dimension n + 1, and is imbedded as a convex 
hypersurface in £ n + 1 . The converse is also true |16| . 
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case there is nothing to adjust, so as for GPEC actions this dropped with increasing A (but was still 
appreciable even for quite large A). 

If we look at the numerical results for Si (S g + XStight) in Table 1 first 
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Table 1 
Results for Si, N = 72 



we can see that Su g ht does indeed drop off with increasing A just like the extrinsic curvature. The behavior 
of S m and max(qi) is also reminiscent of the GPEC action, with the internal geometry becoming more 
regular with increasing A. However, the specific heat C shows only a modest cusp at around A = 2.00 as 
can be seen in Fig. 2, which should be contrasted with the larger peaks seen on these small meshes for both 
the GPEC and Steiner actions. Similarly X2, plotted in Fig. 3, shows no sign of a crumpling transition, 
with only a small increase in the region of the cusp in C, before it rapidly drops off. The value of the 
gaussian term is close to the expected d(N — l)/2, assuring us that the results are equilibrated and both 
the metropolis and flip acceptances are reasonable for all the values of A simulated, so we can be sure that 
the simulation is performing as it should. We also measured the value of the gaussian curvature using 
Sught with the modulus sign removed and found, as expected (our surfaces have spherical topology), 4tt 
for every surface generated. Visual inspection of "snapshots" of the surfaces that arise in the simulation 
confirms the absence of a phase transition, with surfaces looking similar for all of the A values simulated. 
One of these for A = 4.0, but which is typical of all the others, is shown in Fig. 4 and is obviously not 
smooth. Even at the largest A values simulated the surfaces are still some way from satisfying the lower 
bound of 47r on S t i g ht, so it would appear that the disordering effect of S g overcomes Su g ht for all A. 
The behavior of S2 {Area + XSu g ht) is rather bizarre, as can be seen in Table 2. 
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Table 2 
Results for S 2 , N = 72 

Again Su g ht decreases with increasing A but there is no obvious peak in the specific heat. max{qi) now 
increases at large A and there is a huge peak in X2 at around A = 2.5. For A > 2.5 the surfaces generated 
look rather similar to those produced by Si, but near A = 2.5 they are very long, jointed linear structures 
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such as that in Fig. 5. On its own the area action gives surfaces that are collections of long thin spikes 
emanating from a central point, and this is also seen for S2 at small A. It appears that adding Su g ht 
gives (approximately) very long, thin ellipsoids which, while satisfying the convexity property, are not 
the generic smooth surfaces that we envisaged. From the evidence of the specific heat there is little sign 
of a phase transition in this region that might be used to define a continuum limit, though it is always 
possible that a higher order transition may be present. 

Our conclusions are rather disappointing from the point of view of finding candidate random surface 
actions which might be used similarly to the GPEC action to hunt for a non-trivial continuum string 
theory. S\ shows little sign of a transition at all, apart from a modest bump in the specific heat, and S%, 
whilst adding to the bestiary of amusing pathologies that have been observed with various random surface 
models, also looks unpromising. It may of course be possible to incorporate Su g ht as an additional term 
in the GPEC action in order to tune the couplings to see if the approach to the continuum limit could be 
optimized. A further possibility that might be worth pursuing is looking at the effect of self-avoidance, 
which completely changes the behavior of GPEC actions |T^] . 
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Figure Captions 
Fig. 1. The neighbors and angles used in defining Su g ht- 
Fig. 2. The specific heat C for action Si. 
Fig. 3. The gyration radius X2 for action Si. 
Fig. 4. A snapshot of a mesh generated by Si with A = 4.5. 
Fig. 5. A snapshot of a mesh generated by S2 with A = 2.5. 
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